Longitudinal scaling of observables in heavy-ion collision models 
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Longitudinal scaling of pseudorapidity distribution of charged particles (dNch/drj) is observed 
when presented as a function of pseudorapidity (t?) shifted by the beam rapidity {rj - j/bcam) for a 
wide range of collision systems (e"*" + e~ , p+p, d+A and A+A) and beam energies. Such a scaling 
is also observed for the elliptic flow {V2) of charged hadrons in A+A collisions. This is a striking 
observation, as V2 is expected to be sensitive to the initial conditions, the expansion dynamics 
and the degrees of freedom of the system, all of which potentially varies with coUision system and 
colliding energies. We present a study of the longitudinal scalings of dNch/drj, average transverse 
momentum ((px)) and «2 using transport models UrQMD and AMPT for Au+Au collisions at center 
of mass energies (y/s^) of 19.6, 62.4, 200 GeV and Pb+Pb collisions at 2760 GeV. Only the AMPT 
models which includes partonic effects and quark coalescence as a mechanism of hadronization, shows 
longitudinal scaling for dNch/dr], (pt) and 112. Whereas the UrQMD and AMPT default versions 
show longitudinal scaling only for dNch/dr] and (pt). We also discuss the possibility of longitudinal 
scaling of V2 within two extreme scenarios of models with hydrodynamic and coUisionless limits. We 
find the longitudinal scaling of bulk observables to be an important test for the underlying physics 
mechanism in models of particle production. 

PACS numbers: 25.75.Ld 
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I. INTRODUCTION 

Scale invariance in experimental observables from 
heavy-ion collisions at the Relativistic Heavy Ion 
Collider have provided interesting insights about 
particle production mechanism in these reactions [ij . 
Observables like (pt), freeze-out parameters jl], pion 
interferometry radii Q are observed to scale with 
some powers of dNch/drj. Within the frame work 
of thermal models some of these reflect a constant 
energy per particle at freeze-out. The elliptic flow 
for identified baryons and mesons when divided by 
the number of constituent quarks is found to scale 
with the kinetic energy of the particles (or con- 
stituents) This has been interpreted as due to 
development of substantial collectivity in the par- 
tonic phase [5*1 of the evolution of the heavy-ion col- 
lisions and coalescence being the dominant mecha- 
nism of particle production for the intermediate 
range of 2 - 6 GeV/c [g^]. Scalings have been ob- 
served in p-\-p collisions. At low pt (< 2 GeV/c) for 
a given -y/i, the invariant yields of identified hadrons 
are observed t o scale wh en plotted as a function of 
mx (tot = \/p\ — and m is the mass of the 
hadron) , Q • This indicates the absence of significant 
radial flow in p+p collisions. At high pr {> 2 GeV/c) 
the product of the invariant yields of a hadron and 
some power of the ^/s become independent of 
when plotted as a function of (= 2pt/a/s) 
This is interpreted as due to dominance of the pQCD 
process (jets) in p+p collisions. Scaling has been also 
observed in the longitudinal direction, represented 



by a variable rj - yheam, 
directed flow (wi 



V2 [13 and 
ll| in A-l-A collisions. Further 



in dNch/dr] 



such scaling have been widely used to predict the 
values of the observables at hig her energy regimes 
in nucleus-nucleus collisions [12]. Recent results on 
dNch/drj at midrapidity for Pb-|-Pb central collisions 
at y^JNN = 2760 GeV is found to be higher than ex- 
pected from extrapolation based on the longitudinal 
scaling of dNch/drj at lower beam energies [13,]. 

The longitudinal scaling of dNch/dr] is a widely 
discussed subject as it is observed for variety of col- 
liding systems starting from e+ + e~ , p+p, d+K 
to A-l-A collisions [1, [3] ■ This phenomena is of- 
ten called limiting fra gme ntation. It was hypothe- 
sized by Benecke et al. [l^ , Feynman [16|] and Hage- 



dron [17[ that as ^/s — )■ 00 the multiplicity distri- 
bution becomes independent of -^/s- From a mi- 
croscopic picture the longitudinal scaling is under- 
stood, assuming the rapidity distributions of the 
produced particles are functions of x (fraction of the 
hadron longitudinal momentum carried by a typi- 
cal parton) alone, and not of the total energy. If 
the hadron interactions are short-ranged in rapidity 
then the hadron rapidity distributions would repro- 
duce the corresponding distributions of the liberated 
partons. The picture is very similar to Bjorken scal- 
ing of parton distributions. This interpretation can 
be easily linked to initial state gluon dynamics of 
the system [Tsj . From a macroscopic picture, the en- 
tropy conservation in heavy-ion collisions can make 
dNch/dr] insensitive to some aspects of dynamics of 
system and hence may be the cause of the scaling. 
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The pt integrated V2 for a given rapidity range is 
defined as 03 



(1) 



where 4> and 5* are the charged particle azimuthal 
angle and reaction plane angle respectively. The ob- 
served V2 is affected by the initial conditions, it is 
sensitive to the equation of state and its magnitude 
is determined by the interactions of the constituents 
through out the evolution of the system in heavy- ion 
collisions [13, HH ■ Hence the physical interpretation 
of longitudinal scaling of V2 is counter intuitive. Lon- 
gitudinal scaling of V2 exhibited for a wide beam en- 
ergy range, different collision centrality and collision 
species would tend to indicate weak dependence of V2 
on the above mentioned physical conditions. Recent 
studies [211 suggest that the simultaneous observa- 
tion of longitudinal scaling of V2 and dNch/drj can 
be reconciled only if the system formed in heavy-ion 
collisions are weakly coupled. This is contrary to 
other indirect estimations of the shear viscosity to 
entropy ratio which suggests the system is strongly 
coupled [2^ . 

In this paper, we study the longitudinal scaling of 
dNch/dr],{pT) and V2 using models AMPT (A Muhi 
Phase Transport Model, ver 1.11) and UrQMD 
(Ultra Relativistic Quantum Molecular Dynamics, 
ver 2.3) [24] for charged particles in Au4-Au colli- 
sions at y/I^ = 19.6, 62.4, 200 GeV and Pb-|-Pb 
collisions at 2760 GeV. The aim being to see if these 
models also exhibit such longitudinal scalings and 
hence provide a physical insight behind the phe- 
nomena. The UrQMD model is based on a micro- 
scopic transport theory where phase space descrip- 
tion of the reactions and hadron-hadron interactions 
are important. The AMPT model can be studied 
in two configurations, in the AMPT default version 
the minijet partons are made to undergo scattering 
before they are allowed to fragment into hadrons, 
while in the AMPT-SM string melting scenario ad- 
ditional scattering occurs among the quarks and the 
hadronization occurs through the mechanism of par- 
ton coalescence. 



II. LONGITUDINAL SCALING OF dN^i^/dr] 
AND (pt) 




FIG. 1: (Color online) dNch/drj versus ?7-ybeam for 0-6% 
central Au-|-Au collisions at ^snn = 19.6, 62.4 and 200 
GeV from (a) the PHOBOS experiment at RHIC 9], (b) 
UrQMD, (c) AMPT default and (d) AMPT-SM. Also 
shown are the model results from Pb-f-Pb collisions at 
Js^ = 2760 GeV. 
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FIG. 2: (Color online) (pt) versus r;-i/bcam for minimum 
bias Au-f An collisions at ^snn = 19.6, 62.4 and 200 GeV 
from (a) the UrQMD, (b) AMPT default and (c) AMPT- 
SM. Also shown are the model results from Pb+Pb col- 
lisions at v^snn — 2760 GeV. 



Figure [T] shows the dNch/drj versus rj-y\,caia for 0- 
6% central Au-f Au colhsions at ^snn = 19.6, 62.4 
and 200 GeV from (a) the PHOBOS experiment at 
RHIC [9], (b) UrQMD, (c) AMPT and (d) AMPT- 
SM models. Also shown are the results from the 
models for Pb-|-Pb collisions at .^snn — 2760 GeV. 
The ybeam values for ^snn = 19.6, 62.4, 200 and 
2760 GeV are 3.03, 4.19, 5.36 and 7.98 respectively. 
The longitudinal scaling observed in dNc\i / dfj in the 



data (Fig. [Tfa)) is also observed in all the models 
studied. Such scalings have also been observed for 
other models like HIJING ^ and HIJING BB [H, 
27]. 

Figure [5] shows the (px) for the charged parti- 
cles versus ?7-ybcam for minimum bias Au-fAu col- 
lisions at -y/sNN = 19.6, 62.4 and 200 GeV from (a) 
the UrQMD, (b) AMPT and (c) AMPT-SM models. 
Also shown are the results from the models for min- 
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FIG. 3: (Color online) 112 for charged particles versus 
V'yhea.m for 0-40% Central Au+Au collisions at ^snn = 
19.6, 62.4 and 200 GeV from (a) the PHOBOS experi- 
ment at RHIC iia], (b) UrQMD, (c) AMPT defauh and 
(d) AMPT-SM. Also shown are the model results from 
Pb-I-Pb collisions at ^/s^ = 2760 GeV. 



imuni bias Pb+Pb collisions at ^snn = 2760 GeV. 
There are no experimental data available at RHIC 
for (pt) versus rj-yhcum hence not shown in the fig- 
ure. The longitudinal scaling is observed in all the 
models studied. 

These results then sets the stage for studying the 
longitudinal scaling in V2 ■ Note that the goal here is 
not to have a quantitative comparison with data on 
the scalings in dNch/di] and V2, but to see if the ob- 
servations are qualitatively reproduced in the mod- 
els. 



III. LONGITUDINAL SCALING OF V2 

Figure |3] shows the V2 for charged particles ver- 
sus Ty-j/bcam in Au-|-Au collisions at -^snn = 19.6, 
62.4 and 200 GeV The results from the models 
for Pb+Pb collisions at ^snn — 2760 GeV are also 
shown. The collision centrality is 0-40% central and 
is different for that shown for dNch/drj in Fig. [TJ 
The choice of centrality is based on availability of 
the V2 data for charged particles in the experiment 
as a function of rapidity. Figure |3] (a) shows the 
longitudinal scaling of V2 as measured by the PHO- 
BOS experiment 10]. Fig. [3] (b) shows the V2 vs. 
rj-yhciaa from UrQMD model, in (c) the correspond- 
ing results from AMPT default are shown and in (d) 
the same results from AMPT-SM are presented. It 
is observed that the UrQMD and the AMPT default 
models do not show the longitudinal scaling as ob- 
served in the data (Fig. [3] (a)). Only the AMPT 



model with string melting qualitatively reproduces 
the observed longitudinal scaling of V2. 

IV. DISCUSSIONS 

It is worthwhile to now discuss briefiy the differ- 
ences in these transport models. The main differ- 
ence between UrQMD and AMPT lies in the initial 
conditions (for AMPT taken from HIJING (H) and 
additional initial state rescatterings in AMPT. The 
main difference between AMPT default and AMPT- 
SM lies in the following: The string melting ver- 
sion of the AMPT model is formulated on the idea 
that for energy densities beyond a critical value of 
~ 1 GeV / fm^ , it is difficult to visualize the coex- 
istence of strings (or hadrons) and partons. Hence 
the need to melt the strings to partons. This is done 
by converting the mesons to a quark and anti-quark 
pair, baryons to three quarks etc. The scattering 
of the quarks are then carried out based on par- 
ton cascade [1^. The parton-parton cross section 
taken here is 10 mb. Once the interactions stop, 
the partons then hadronizes through the mechanism 
of partonic coalescence. While for the AMPT de- 
fault the scattering occurs for minijet partons (no 
melting of strings to partons) and hadronization oc- 
curs through fragmentation process [2^ . This model 
based study then suggests that partonic interactions 
in high energy density matter is essential to quali- 
tatively reproduced the simultaneous observation of 
the longitudinal scalings in dNch/di] and V2 in ex- 
periment. If this is the actual cause then it will be 
interesting to have experimental measurements of V2 
vs. 7] for lower beam energies where we do not expect 
to create a sufficiently high energy density system to 
see the breakdown of such a V2 longitudinal scaling. 

We now briefly discuss some other possibilities 
which could explain the longitudinal scaling of V2. 
One of them is based on the arguments whether the 
system is weakly coupled or strongly coupled. A 
weakly coupled system has been argued to favor the 
combined V2 and dNch/drj scaling behavior [2l|. It 
has been suggested that for systems where the inter- 
actions among the constituent particles are small, or 
the system is close to free streaming, called the col- 
hsionless limit HI, the ?;2 ~ ^ Vr^- Where v 
is the relative velocity of the particles, a is the mo- 
mentum transfer interaction cross section and the 
product nRxRy is the transverse overlap area for 
the two nuclei. In this model, one can easily see 
that V2 should exhibit a longitudinal scaling simi- 
lar to dNch/dr] provided (va) does not change with 
beam energy. A linear dependence of V2 with change 
in jj;^ ^ has been observed in experiments over 
a wide collision systems [3(i]- In the event of (va) 
changing with beam energy, possibly due to change 
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in the relevant degrees of freedom (hadronic or par- 
tonic), there would be a breakdown! of the longi- 
tudinal scaling of V2- This is consistent with the 
conclusions from our model study. Now let us move 
to the other extreme limit, where the rescattering 
among the consitutent particles are abundant lead- 
ing to the hydrodynamic limit [29j . In such a model 
the V2 is proportional to the average transverse mo- 
mentum of the particles among several other quan- 
tities as discussed in [l^. If the (pt) also exhibits 
a longitudinal scaling then V2 in the hydrodynamic 
limit scenario should also exhibit the scaling. Mea- 
suring {pt) vs. T] could help address the cause of the 
longitudinal scaling of V2. However we have seen in 
Fig. [5] that the models based on transport approach 
also exhibit longitudinal scaling of (pt)- The model 
study for all the three observables indicates that ob- 
serving longitudinal scaling in dNch/drj and/or (px) 
does not neccessarily implies we should see a similar 
scaling in V2- 

V. SUMMARY 

In summary, we have discussed the simultaneous 
observation of longitudinal scaling of dNch/drj and 
V2 when plotted as a function of jy-ybeam in RHIC 
experiments in Au-|-Au collisions at ^snn = 19.6, 
62.4, and 200 GeV. There are several physical argu- 
ments based on convervation of entropy and hadron 
distributions in rapidity corresponding to the rapid- 
ity distributions of partons (assuming hadron inter- 
actions are short ranged in rapidity) to explain the 
longitudinal scaling of dN^h/dii. The UrQMD and 
AMPT models qualitatively reproduce the the ex- 
perimental observations for the .^snn — 19.6 - 2760 
GeV. These models also exhibit longitudinal scal- 
ing in (pt)- The observation of longitudinal scal- 
ing for V2 which could be senstive to several quani- 
ties like initial conditions/densities, equation of state 
and rescatterings among the consituents is bit in- 



triguing. We find that among the models studied, 
only AMPT with string melting qualitatively shows 
the behaviour as exhibited by the data. The main 
difference between this version of the AMPT model 
compared to the default version and UrQMD lies in 
the partonic interactions and hadronization through 
coalescence mechanism. We also discussed that the 
longitudinal scaling can naturally occur for a system 
in the collisionless limit. In such a limit V2 propor- 
tional only to dNch/dr] provided the product of inter- 
action cross section and average relative velocity of 
the particles does not change with the beam energy 
studied. It is expected that the interaction cross 
sections could be very different if the relevant de- 
grees of freedom are partonic or hadronic. So study- 
ing the logitudinal scaling of V2 at lower and higher 
beam energies will provide additional insight. For 
the other extreme, hydrodynamic limit, since V2 is 
proprotional to (pt), it would be interesting to mea- 
sure in experiments the (px) vs. rj for various beam 
energies. However the model study of the obsrvables 
dNch/drj, (pt) and V2 indicates that observing lon- 
gitudinal scaling in dNch/drj and/or (px) does not 
neccessarily implies that such a scaling will follow in 
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